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ABSTRACT 

PSR J1833-1034 and its associated Pulsar Wind Nebula (PWN) has been investigated in depth 
through X-ray observations ranging from 0. 1 to 200 ke V. The low energy X-ray data from 
Chandra reveal a complex morphology that is characterised by a bright central plerion, no 
thermal shell and an extended diffuse halo. The spectral emission from the central plerion 
softens with radial distance from the pulsar, with the spectral index ranging from r = 1.61 in 
the central region to T =2.36 at the edge of the PWN. At higher energy INTEGRAL detected 
the source in the 17-200 keV range. The data analysis clearly shows that the main contribution 
to the spectral emission in the hard X-ray energy range is originated from the PWN, while the 
pulsar is dominant above 200 keV. Recent HESS observations in the high energy gamma-ray 
domain show that PSR J1833-1034 is a bright TeV emitter, with a flux corresponding to ~2 
per cent of the Crab in 1-10 TeV range. In addition the spectral shape in the TeV energy region 
matches well with that in the hard X-rays observed by INTEGRAL. Based on these findings, 
we conclude that the emission from the pulsar and its associated PWN can be described in a 
scenario where hard X-rays are produced through synchrotron light of electrons with Lorentz 
factor y ~10 9 in a magnetic field of ~10 micro Gauss. In this hypothesis the TeV emission is 
due to Inverse Compton interaction of the cooled electrons off the Cosmic Microwave Back- 
ground photons. Search for PSR J1833-1034 X-ray pulsed emission, via RXTE and Swift 
X-ray observations, resulted in an upper limit that is about 50 per cent. 
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rays: observations 



1 INTRODUCTION 

G2 1.5-0.9 is a Crab nebula like supernova remnant (SNR) with a 
plerionic structure, i.e. with a filled center form and a flat radio 
spectrum F y oc V with a between and -0.3 (Weiler & Panagia 
1978). It does not display any clear evidence of a thermal shell. 
Through radio observations, the system was found to include a 61.8 
ms pulsar near to the centre of the nebula (Camilo et al. 2006). The 
X-ray morphology investigated in depth through several Chandra 
observations (Matheson and Safi-Harb, 2005, Safi-Harb et al. 2001, 
Camilo et al. 2006) shows a near spherical structure with a bright 
compact source surrounded by the PWN (~40 arcsec in radius), 
which is in turn surrounded by non- thermal halo of about 150 arc- 
sec in radius. This extended halo shows evidence of a possible SNR 
structure, exhibiting filamentary-like features and limb brightening. 
However, the X-ray spectral emission from this region is decidedly 
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non thermal in character, as are the brighter knots within its en- 
velope. The origin of this component is still unclear, but its non- 
thermal nature suggests it is an extension of the synchrotron nebula 
(Safi-harb et al. 2001), with a contribution of dust scattered X-rays 
from the plerion (Bocchino et al. 2005). The detailed X-ray Chan- 
dra image of the central zone shows that the compact source, i. 
e. the site of the pulsar, is surrounded by a double lobed ellipti- 
cal emission region of size ~7 arcsec by 5 arcsec (Camilo et al. 
2006). The offset of the pulsar from the geometrical centre of the 
G2 1.5-0.9 shell is <5 arcsec. The elliptical geometry is possibly an 
indicator of a toroidal structure around the pulsar as observed for 
the Crab and Vela pulsars or eventually the putative axisymmet- 
ric jet structures observed by Fiirst et al. (1988) by means of high 
resolution 22.3-GHz observations of this system. 

The X-ray emission, as observed by Chandra, is well de- 
scribed by a power law model with a spectral index that is a func- 
tion of the angular distance from the central region (Safi-Harb et 
al. 2001). The point-like X-ray source within a 1 arcsec region at 
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the site of the pulsar has a spectral index of T = 1.4. The spectral 
emission from the central plerion softens with the radial distance 
from the pulsar, and shows a spectral index ranging from a value 
T = 1.61, in a region within the central 5 arcsec, to T = 2.36, at 
the edge. In the 0.1-10 keV energy range the outer halo emission 
is generally well described by a softer spectrum with T ~2.4, while 
the knots have spectral indices in the range from T ~2. 14 to T ~2.5. 
In 0.5 - 10 keV energy range, the total unabsorbed flux from the 
extended component and compact core is 1.1 XlO -10 erg cm 2 s -1 
and the corresponding N H value is 2.24 xlO 22 cm 2 . Following the 
detailed discussion presented in Camilo et al. (2006) we adopt a 
distance estimate of 4.7 ± 0.4 kpc. 

Careful searches in the Chandra and XMM data has failed to 
reveal any pulsed X-ray emission at the period and ephemeris of 
the radio pulses. La Palombara and Mereghetti (2002) set an upper 
limit on the 0.5 to 10 keV pulsed fraction between 7.5 per cent and 
40 per cent, however, they did not search pulsed events in the range 
of frequencies that comprises the spin frequency predicted by the 
radio ephemeris (the pulsed emission was discovered three years 
later by Camilo and co-workers). With a period of 61.8 ms and a 
period derivative of 2.02 x 10~ 13 s s _1 , the characteristic spin down 
age for a braking index of n = 3 is 4.8 ky. PWN evolution in the 
light of the spinning neutron star energy loss rate, suggests the age 
should be much less that 5 ky (Chevalier 2004). There is good his- 
torical evidence from ancient Chinese records to believe that PSR 
J1833-1034 is associated with a guest star supernova explosion that 
took place in BC 48, making the system just over 2050 years old 
(Wang et al., 2006). This younger age implies the pulsar was born 
with a period close to 47 ms. 

Finally, very recently HESS telescope has detected Very High 
Energy (VHE) emission in the TeV energy range from PSR J1833- 
1034 and its associated PWN (Djannati-Atai et al. 2007). 

In this paper we present the analysis of the INTEGRAL-IBIS 
(Ubertini et al. 2003) hard X-ray observation of PSR J1833-1034. 
The spectral information about this pulsar and its associated PWN 
above 10 keV was up to now still missing. To search for X-ray 
pulsation we performed a detailed temporal analysis of the avail- 
able archival RXTE and Swift-XRT data. In Sect. [2] we inves- 
tigate the hard X-ray energy spectrum of the pulsar and PWN 
through INTEGRAL observation, while the TeV HESS emission 
from PSR J 1833- 1034 and its association with the hard X-ray coun- 
terpart is discussed in Sect. [3] Timing analysis is presented in Sect. 
|4]and finally our conclusions are drawn in Sect. [5] 



2 THE HARD X-RAY SPECTRA OF THE PULSAR AND 
PWN OBSERVED BY INTEGRAL 

PSR J1833-1034 has been observed as part of the third INTEGRAL 
IBIS/ISGRI (here after INTEGRAL) survey processing (Bird et al. 
2007), which consists of all exposures from the beginning of the 
mission (November 2002) up to April 2006. The total exposure 
on PSR J1833-1034 is of ~ 1.7 Ms. INTEGRAL images for each 
available pointing were generated in various energy bands using 
the ISDC offline scientific analysis software OSA (Goldwurm et al. 
2003) version 5.2. The individual images were then combined to 
produce a mosaic of the entire sky to enhance the detection signif- 
icance using the system described in detail by Bird et al. (2007). 
Fig. [U shows the 17^-0 keV energy band image of the region sur- 
rounding PSR J 1833- 1034 and a clear excess is detected at a signif- 
icance of ~25<x. The INTEGRAL image is consistent with a point- 
like source that is coincident with the site of PSR J1833-1034. The 



Figure 1. INTEGRAL image of PSR J1833-1034 in 17^10 keV range. The 
spatial extension of the excess reveals a point-like source. The circle repre- 
sents the extension of HESS J 1833- 1034 (5 arcmin radius). In the insert we 
show the INTEGRAL spectrum. 



data are well fitted by a simple power-law with photon index T = 
2.2 +0 \ (the spectrum is shown in the insert in Fig. QJ. The 20-100 
keV flux of 5.2x10"" erg cm" 2 s -1 , as measured by INTEGRAL, 
comprises 0.4 per cent of the pulsar spin down luminosity of E = 
3.38 xlO 37 erg s~', assuming a distance to the source of 4.7 kpc. 
The most likely site of this emission should be from the central 
bright plerion. Chandra observations in 0.1-10 keV energy range 
have clearly detected the softening of the X-ray emission with ra- 
dial distance from the pulsar, that is presumably the source of ener- 
getic electrons. This evidence has all the hallmarks of synchrotron 
cooling, so that we would expect the hard X-rays as observed by 
INTEGRAL to originate nearby the central source. 

To evaluate the contribution to the INTEGRAL spectrum of 
the different components of the source (pulsar, PWN and extended 
emission), we plot in Fig.|2]for each component the spectral model 
as observed by Chandra (Safi-Harb et al. 2001). The low surface 
brightness extended emission (originated between 50 arcsec and 
150 arcsec radius around the pulsar), is reproduced with a soft 
power-law (r=2.36, dashed line). The emission of the pulsar (dot- 
dashed curve), extracted from the central innermost region in a 
circle with 1 arcsec radius, is reproduced with a hard power-law 
(T=1.4). The emission of the PWN (dotted line) is extracted from 
a circle of 40 arcsec radius and is best fitted by a softer power-law 
with respect to the pulsar (r=2.3). All components are absorbed 
by a cold gas with column density N H =2.2xl0 22 cirT 2 as derived 
through Chandra (Safi-harb et al. 2001) and XMM-Newton (War- 
wick et al. 2001) observations. In Fig.[2]fhe INTEGRAL data are 
also shown. From this figure it appears the pulsar contributes very 
little in the 20-200 keV energy range, and the extended emission 
is too weak with respect to the other components to affect the high 
energy emission as detected by INTEGRAL. Therefore the PWN is 
likely to contribute most to the hard X-ray flux. Indeed, the contri- 
bution of the pulsar at the 20-100 keV flux, as derived through IN- 
TEGRAL data, is lower than 20 per cent and it become even lower 
in the soft X-ray energy range, i.e. ~ 3 per cent (see also Kargaltsev 
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Figure 2. Spectral Energy Distribution from soft X-ray to TeV energy 
range. Each spectral component identified through Chandra observation 
has been plotted separately. INTEGRAL data are plotted on the top of the 
model. We also plot the Fermi-GST (formerly GLAST) LAT 1-year sensi- 
tivity curve taking into account the diffuse Galactic and the residual instru- 
mental background. 

& Pavlov 2008). This is clearly visible in the SED plotted in Fig.[2] 
In 20-100 keV energy range a fraction of the total emission that is 
larger of 80 per cent is originated in the PWN. Using the pulsar spin 
down luminosity of E = 3.38 x 10 37 erg s _1 , the total energy conver- 
sion efficiency (in 20-100 keV energy range) 
per cent, with the PWN contribution that is L PWN /E ~ 0.3 per cent 
and the pulsar L puhar jE ~0.1 per cent. The value of the total en- 
ergy conversion efficiency r\ m is well in agreement with the aver- 
age value measured in the others pulsar/PWM systems observed by 
INTEGRAL, that is around 1 per cent (PSR J 1846-0258, McBride 
et al., 2008; PSR J1617-5055, Landi et al., 2007; Vela, Hoffmann 
etal., 2007; PSR J1513-5906, Forot et al., 2006; PSR J181 1-1925, 
Dean et al., 2008). PSR B0540-69 (Campana et al. 2008) is char- 
acterized by the highest value rj m ~ 5.8 per cent, while Vela by 
the lowest r\, a ~ 0.02 per cent. Finally we note that the ratio be- 
tween the energy conversion efficiency of the PWN and pulsar is 
T) pwn i [ T) psr ~ 4-3, in very good agreement with the mean value (~4) 
for a sample of X-ray pulsar/PWN systems observed in the soft X- 
ray energy range with Chandra (Kargaltsev & Pavlov 2008). 



3 ASSOCIATION WITH HESS J1833- 105 

Observations in the TeV energy range have demonstrated that 
PSR J1833-1034 and its associated PWN are powerful emitters in 
the VHE gamma-rays (Djannati-Atai et al. 2007). The position of 
HESS J1833-105 (18*33 m 32.5 s ± 0.9 s , -10d33'19" ± 55") is in very 
good agreement with the location of the source detected by IN- 
TEGRAL. The properties of the hard X/gamma-ray emission from 
the pulsar/PWN system observed by INTEGRAL and HESS are re- 
ported in Table [T] while in Fig. Q] the HESS position is shown as 



a circle on the INTEGRAL image. The photon index of the TeV 
spectrum is 2.1 ±0.2 with a flux corresponding to ~2 per cent of 
the Crab (see Table[TJ. In Fig. [2] we report the best-fit power-law 
reproducing the spectrum of HESS J1833-105, with the three lines 
representing the upper, lower and best fit value of the slope. The 
shape of this emission, extrapolated in the X-ray energy range, 
matches well with the PWN spectrum observed by INTEGRAL, 
with the ratio between the X-ray and gamma-ray luminosity of 
L x /L r ~ 30. We assume that the soft-gamma ray emission is pro- 
duced through synchrotron process of electrons in a magnetic field 
of 10 micro Gauss. The value of the magnetic field is close to that 
implied from gamma-ray observations assuming (a) the gamma- 
ray flux is produced by Inverse Compton (IC) scattering on Cosmic 
Microwave Background (CMB) photons and (b) a ratio h x fL y ~ 
30 (Djannati-Atai et al. 2007). The equipartition value suggests 
the presence of a magnetic field much higher, of about B=0.3 mG 
(Camilo et al. 2006). However, HESS observations seems to sug- 
gests that the magnetic field in PWN can be much lower than the 
equipartition value (Chevalier 2004). To produce the observed 20- 
100 keV emission through synchrotron light in 10 micro Gauss 
magnetic field, we need electrons with Lorentz factor of ~10 9 . 
This value remains still very high even if we assume a higher mag- 
netic field, e.g. y ~10 8 for B=50 micro Gauss. In this scenario TeV 
emission can be produced through IC scattering of the same elec- 
trons population with CMB light. For electrons with Lorentz fac- 
tor of ~10 9 , the maximum energy gained in a single scattering is 
v max = y 2v CMB ~ 10 14 eV. The energy of the electrons producing 
TeV emission can be lower with respect to that required to produce 
the 20-100 keV spectrum. This matches well with the hypothesis 
that the TeV emission is produced by the electrons cooled through 
synchrotron radiation in a 10 micro Gauss field. 

Future Fermi-GST observations will allow a detailed study of 
the SED in the MeV-GeV energy range, where the pulsar contribu- 
tion should be dominant. The 1 year LAT sensitivity curve is also 
plotted in Fig. [2] 

In the sample of pulsar/PWN systems detected by INTEGRAL, 
PSR J1833-1034 together with Crab and PSR J1846-0258, is char- 
acterised by a high Lio-ioofev/Li-ioreV ratio and by positional co- 
incidence with HESS. This high ratio seems to be an artefact of 
young pulsar/PWN systems for which the TeV source is coincident 
with the pulsar, and maybe due to the fact that the electrons (syn- 
chrotron cooled) reservoir is not completely filled because of the 
lack of accumulation time. 



4 TIMING ANALYSIS 

In the case of INTEGRAL data, the combination of the cumula- 
tive long observation period and the relatively low statistical sig- 
nificance of the detection make impossible to isolate and measure 
a pulsed hard X-ray component in the emitted flux. We then per- 
formed a period search using RXTE/PCA data in order to estab- 
lish an upper limit on the pulsed emission in the range 2-10 keV. 
PSR J1833-1034 has been observed once by RXTE, during pro- 
posal ID 20259 performed in November 8, 1997 with G21.5-0.9 as 
a target. We have selected PCA Good-Xenon observations, i.e. the 
data mode with the full timing resolution. The data were filtered 
using the standard criteria (elevation angle, south atlantic anomaly 
passages, etc); the baricentrisation was performed using the task 
f xbary and the pulsar location from Chandra (Camilo et al. 2006). 
The total exposure time on the source was 73 ks. We performed a 
period search in a reasonably large frequency range around the ex- 
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Table 1. Hard X/gamma-ray properties of PSR J1833-1034 and its PWN detected by INTEGRAL and HESS 
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trapolation of the radio ephemeris to the RXTE observation time, 
since the radio ephemeris were not available at the time of RXTE 
observations. We used the Z, 2 test (Buccheri et al. 1993) with m = 2 
and m = 10 harmonics. We ignored the Z 2 test with m = 1 because 
this is very powerful for sinusoidal peaks, while we have no idea 
on lightcurve expected in our case (Protheroe 1987). We found no 
pulsation at the 99 per cent confidence level. Upper limits on the 
pulsed emission were derived with the hypothesis that the pulse 
shape is either sinusoidal or very narrow (in the limit of a Dirac 
delta shape), following the prescriptions of Protheroe (1987) and 
De Jager et al. (1989). For the 2-10 keV band we obtained a flux 
value between 2.5 and 1.2xl0~ 5 photons cirr 2 s~' for a sinusoid 
and a Dirac delta pulse, respectively. This pulsed flux upper limit 
corresponds to 1.7x 10~ 13 erg cirT 2 s~' and 8.8xl0~ 14 erg cirT 2 s~' 
respectively, assuming that the pulsar photon index is T = 1.5 as 
found by Chandra. These upper limits are of the same order of 
magnitude, if not slight more stringent, than those set by Chandra 
(Camilo et al. 2006). The pulsed emission is therefore lower than 
50 per cent of the total emission from the pulsar. 

We also analysed observations of PSR J1833-1034 made with 
the XRT telescope onboard the Swift satellite (Gehrels et al. 2004). 
We selected all the observations having this source in the field of 
view and performed in the Windowed Timing (WT) mode, that en- 
sures the full timing capability of this instrument, with a time reso- 
lution of about 1.7 ms. The dataset contains 9 pointings, performed 
between October 2 and 21, 2007. The total net exposure time af- 
ter the standard pipeline filtering is about 52 ks. We performed a 
period search using an approach similar to the one used for RXTE 
observation, and still no significant pulsation was found at the 99 
per cent confidence level. The upper limit on the pulsed flux was of 
~5xlO~ 13 erg cnT 2 s~'. 



5 CONCLUSIONS 

The combined spectral study of the Crab-like pulsar PSR J1833- 
1034 and its associated PWN over the Chandra and INTEGRAL 
energy range, strongly suggests that the 20-100 keV flux is mainly 
due to the PWN component (more than 80 per cent of the total 
flux), and the PSR contribution can only be dominant above 200 
keV. PSR J1833-1034 is also a bright TeV emitter as detected by 
HESS. The INTEGRAL emission represents the hard X-ray coun- 
terpart to the VHE HESS source and the TeV emission matches 
well with the 17-200 keV spectral shape. But for PSR J1617- 
5055 and PSR J15 13-5906, most of the pulsar/PWN systems ob- 
served by INTEGRAL have the hard X-ray emission dominated by 
the PWN. However, a detailed comparison between these sources 
shows that the relative contribution from the pulsar and the PWN 
differs very much from source to source. The observed properties 
of PSR J1833-1034 suggest the hard X-ray emission is produced 
through synchrotron light of electrons with y ~10 9 in a magnetic 
field of ~10 micro Gauss while the TeV emission is due to IC of 
the cooled electrons on the CMB photons. 

We note that in the sample of pulsar/PWN systems observed 



by INTEGRAL PSR J1833-1034 together with Crab and IGR 
J1846-0258 are the only systems showing similar properties, i.e. 
high L 2 o-iookev/Li-ioTeV ratio (»1 ) and with the position of the 
HESS coincident with that of the pulsar/PWN. This evidence sug- 
gests that the cooled electrons responsible for the TeV emission 
remain confined in the same region of the synchrotron emission. 
This could be an artefact of the young pulsar/PWN system. 

We did not detected any evidence of pulsations in RXTE and 
Swift-XRT data (with upper limit for the pulsed emission of 50 per 
cent of the total). The lack of X-ray pulsed emission is not fully 
unexpected in view of the moderate magnetic field of the pulsar, 
B=3.58xl0 12 Gauss. However, we stress also that Crab and PSR 
B0540-69, clearly exhibiting X-ray pulsations, are characterized by 
a magnetic fields of 3.78xl0 12 and 4.96xl0 12 Gauss respectively, 
i.e. values similar to that observed in PSR J1833-1034. The lack 
of X-ray pulsed emission from PSR J1833-1034 could also be 
explained with the radio and X-ray beams pointing to different 
directions or with a different aperture, and the line of sight sweeps 
only one of them. A similar hypothesis has been proposed to 
explain the behaviour of Geminga and the radio-quiet pulsars 
(Harding et al. 2007). 

This program is funded by Italian Space Agency grant via con- 
tracts 1/008/07/0 and 1/088/06/0. 
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